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ABSTRACT 

We have examined the recently detected very high energy (VHE) pulsed radiation from the 
Crab pulsar. According to the observational evidence, the observed emission (> 25GeV) peaks at 
the same phase with the optical spectrum. Considering the cyclotron instability, we show that the 
pitch angle becomes non-vanishing leading to the efficient synchrotron mechanism near the light 
cylinder surface. The corresponding spectral index of the emission equals —1/2. By studying 
the inverse Compton scattering and the curvature radiation, it is argued that the aforementioned 
mechanisms do not contribute to the VHE radiation detected by MAGIC. 



Subject headings: 
non-thermal 

Introduction 



instabilities - plasmas - pulsars: individual (PSR B0531+21) - radiation mechanisms: 



For studying the high energy radiation of 
pulsars, two major m echanisms are used: the 



synch rotron process (Pacini 1971 : Shklovskvl 



1970 ) and the Inverse Compton Scattering (ICS) 



(Blandford ct al. 1990). lowever, the problem 
of identification of a location, where the radiation 
comes from still is a matter of discussion. For solv- 
ing this problem an appr oach of the so-ca lled polar 
cap model was proposed (|Sturrocklll97l[ h Accord- 
ing to this model, due to very strong electrostatic 
field, particles are uprooted from the star's sur- 
face l ayer and accelerate along the m agnetic field 
lines ( Ruderman fe Sutherland 19751 ) inside a zone 
where the electric field is nonzero (gap). This 
process leads to relativistic energies of electrons, 
which in turn, create radiation. Unfortunately, the 
particle energy inventory accumulated within the 
gap is not enough to explain the observed high en- 
ergy emission. For solving this problem, to enlarge 
somehow the gap zone, series of works have been 



invoked. Arons fe Scharlemann (1979) have con 



sidered an effect of rectifying the curved magnetic 
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field lines and they have shown th at this process 
l eads to an increase of the gap size. Harding et alj 
(198 l|) have applied t he me thod developed by 
Arons fe Scharlemann ( 19791 ) for studying the 



high altitude emission from the pulsar slot gaps. 
A three-dimensional model of optical and 7-ray 
emission from the slot gap accelerator was exam- 
ined and it has been predicted that the emission 
below 200MeV must exhibit correlations in time 
and phase with the radio band. Somewhat differ- 
ent mechanism leading to the effect of the gap size 
increase was proposed by lUsov fc Shabadl (| 19851) 
where the intermediate formation of positronium 
(electron-positron bound state) was studied. Gen- 
erally speaking, the space-time close to neutron 
star' s surface is slightly c urved . Based on this 
fact, iMuslimov fc Tsvganl (| 19921 ) have considered 
an effect of creation od an additional electric field 
due to the influence of the Kerr metric. All afore- 
mentioned mechanisms cannot provide a signifi- 
cant increase of the gap size to explain the ob- 
served high energy radiation from pulsars. To 
solve this problem, the so-called out er gap model 



has be en proposed [see for example (jCheng et al 
1986al lbl)]. According to this approach radiation is 



formed further out in the magnetosphere of pul- 
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sars and mechanisms responsible for emission are: 
the synchrotron, inverse Compton and curvature 
radiation processes respectively. 

Recently, MAGIC Cherenkov telescope has de- 
tected from the Crab pulsar the VHE pulsed 
emission above 25GeV (lAliu et all |2008( ). The 



most peculiar feature of the observed radiation 
is coinc idence of t he h igh energy and optical 
signals (jAliu et al.l 120081 ). This indicates that 



the polar cap models cannot be applied for ex- 
plaining the MAGIC observ ational data. In 
(|Machabeli k Osmanovll2009l ) we have examined 



the Crab pulsar's recently observed VHE pulsed 
emission data. It is worth noting that due to very 
small cooling timescales, particles rapidly transit 
to the ground Landau state completely killing the 
subsequent radiation. We have found that due to 
the cyclotron instability, the optical spectrum is 
generated, which in turn via the quasi-linear dif- 
fusion provokes the increase of the pitch angles, 
leading to the synchrotron process with the spec- 
tral index —1/2. It has been shown that the emis- 
sion in the optical and high energy (> 25GeV) 
bands originate from the well localized regions, 
leading to the observational fact that the signals 
peak with the same phases. As we have already 
noted, the analysis of the observational data indi- 
cates that the curvature radiation and the inverse 
Compton mechanism must be excluded from the 
possible emission mechanisms. Generally speak- 
ing, for very strong corotating magnetic fields, the 
particles a) move along curved trajectories and b) 
accelerate up to relativistic energies. Therefore, 
it is of great importance to understand why the 
curvature radiation and the ICS are not involved 
in the process of emission. 

The paper is organized as follows. In Section 
2 we consider the synchrotron, Inverse Compton 
and curvature radiation mechanisms, in Sect. 3 we 
present our results and in Sect. 4 we summarize 
them. 

2. Main consideration 

In this section we are gong to consider three ma- 
jor emission mechanisms: the synchrotron process, 
the ICS and the curvature radiation respectively. 



2.1. Synchrotron emission 

As we have already mentioned, timescales of 
transit of particles to the ground Landau state 
is very small and electrons very soon start mov- 
ing al ong the magnetic field lin es without emis- 



Machabeli fc Usovl (|l979h have found that 



the situation changes due to the cyclotron instabil- 
ity, leading to creation of nonzero pitch angles and 
the correspon ding radiation process. Th e method 
developed in dMachabeli k Usovl (119791)) was ap- 
plied by us in (|Machabeli k Osmanovil2009l ). We 
have shown that the quasi-linear duffusion excites 
the transverse and longitudinal-transverse waves: 



u t = kc(l- Si) , 
wit = k c(l - Si - 5 2 ) 



(1) 
(2) 



where k is the modulus of the wave vector, k.. and 
k± are the wave vector's longitudinal (parallel to 
the magnetic field) and transverse (perpendicular 
to the magnetic field) components respectively, c 
is the speed of light and 



Si = 



,2 Jl 



(3) 



were, u p = \J Aim p e 2 j 'm is the plasma frequency, 
lub = eB/mc is the cyclotron frequency, e and 
m are electron's charge and the rest mass respec- 
tively and n p is the plas ma density. According 



to a m odel developed by iMachabeli fc Osmanov 
(2009), we consider the plasma composed of two 
components: a) the plasma component with the 
Lorentz factor, 7 P and b) the beam component 
with the Lorentz factor, 7f>. 

( 19921) 



Kaz begi et al 



have shown that the 
aforementioned modes generate if the cyclotron 
resonance regime 



&ll k x ii x i — — — 0, 



(4) 



takes place. u x = cV jb/ P^>b is the drift velocity, 
V.. is the longitudinal velocity component and p is 
the curvature radius of field lines. 

For the parameters of the Crab pulsar, P « 
0.033sec R s w 10 6 cm, n ps « 1.4 x 10 19 cm" 3 , 
B s « 7 x 10 12 G, 7 b w 10 8 , from Eqs. (H) one 
can show that the development of the cyclotron in- 
stability occurs in the optical band (~ 10 15 ) close 
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to the light cylinder zone. P denotes the pulsar's 
period, R s - its radius, and n ps and B s are the 
plasma density and the magnetic field induction 
respectively, close to the star. 

When electrons emit in the synchrotron regime, 
they experience the so-called radiative forces 



(|Landau fc Lifshitzlll97lh : 



F\ 



-o^(l + 7 b > ) 



2 K\ 



where ip is the pitch angle. These forces try to 
decrease the pitch angle contrary to the quasi- 
linear diffusion, which tends to widen the value 
of ip. The dynamical process saturates when 
these two processe s balance each other. As it 
has been shown (Machabeli fc Osmanov .2009; 
Malov fc Machabelil l200l[) . when 'yip » 1 the 
kinetic equation in the quasi-stationary regime 
(d/dt — 0) reduces to 



d[Fj] i d[$F ± f\ _ 1 d 2 



p..ip dip 



ip dipdp,. 



ld_ 

ip dip 



f a a 

\ dip ~ L " dp 



n 21 

x » dip 



f 



(6) 



where / = f(ip,p\\) is the distribution function of 
particles, p.. is the longitudinal momentum, 



D 



9 9 

it e ribC 
2lu ' 



D 



ir 2 e 2 nbUJB 
2mcy 2 uj 2 



. (7) 



are the diffusion coefficients and rib = B/(Pce) 
is the density of the beam component. This 
equation can be solved easily if one expresses 
the distribution function as x(V0/(Pii )• Then, 
after substituting t his expression into E 



Chkheidze fc Machabelill2007l ): 
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one gets following (Machab eli fc Osmanovl 1200' 



xW = C 1 e-^,/(p || ) = 



C* 2 



(/ 9 9 T7 2 e 2 ibnhC 
a w % 



lb 



where 



and 



A 



4e 6 B 4p3 7 4 

37r 3 m 5 c 7 76 ' 



(8) 
(9) 



r °° iPxWdip _ 0.5 



Trr- ( 10 ) 



is the mean value of the pitch angle. 

Relativistic electrons moving in the magnetic 
field emit photons with en ergies expressed by e ~ 
1.2x lQ- 17 ff7 2 sin^(Gey) (|e.g. Rvbicki fc Lightman 
1979), which after applying Eq. (fTOjl reduces to: 



e syn (GeV) « 6 x 10" 



18 f3^m^jl 
4P 3 e 6 7 4 



(11) 



Accord ing to the results of the MAGIC Cherenkov 
telescope (|Aliu et al.ll2008h . the VHE (> 25GeV) 
pulsed radiation from the Crab pulsar peaks with 
the same phase as the optical signal. This circum- 
stance is in a good agreement with our model, 
since, as we have shown, the pitch angles increase 
due to the cyclotron instability, which in its turn 
occurs in the optical frequency ranges. 

Let us apply Eqs. (|10ll 1|) to the Crab pulsar. 
Then by taking into account the following param- 
eters, R s w 10 6 cm, B s fa 7 x 10 12 G and ~f p ps 3, 
one can show that the pitch angle is of the order 
of 10 -5 , which guarantees the observed high en- 
ergy emission if 7& ps 3.2 x 10 8 . For the mentioned 
set of parameters one has a^ 2 7 2 S> ir 2 e 2 ipnbc/jb, 
which reduces the distribution function to /(p ) cx 



On the other hand, the spectrum of the 
synchrotron rad iation behaves as I v cx v~^~ 
(jGinzburgl Il981l ). where (5 describes the particle 
distribution function, / cx 7^. This indicates 
that, for the case considered in the present paper 
((3 = 2) the spectral index of the VHE synchrotron 
emission equals —1/2. 

Therefore, as we see, the synchrotron emission 
can explain the observed VHE radiation, and as 
we have already seen, for this purpose the par- 
ticles must have very high Lorentz factors. On 
the other hand, these particles will inevitably en- 
counter soft photons, which in turn can also create 
the high energy radiation via the ICS. But in this 
case the emission will not be localized contrary 
to the observational evidence, indicating that for 
some reason the ICS is not involved in the process 
of the detected emission. The next subsection is 
dedicated to this particular problem. 

2.2. Compton scattering 

It is well known that when a photon with energy 
e encounters a relativistic electron, under certain 
conditions photons might gain energy. The cor- 
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responding frequency after scat tering is given by 
(|e.g. Rvbicki fc Lightman1ll979l ): 



l-/3cos6> 



l-/3cos(9' + -^(l-cos6>' 



(12) 



where lo is the frequency before scattering, (3 = 
v/c, 6 = (PK), 0' = (PK'); 6" = (KK')- % K 
and K' we denote the three momentum of the pho- 
ton before and after scattering respectively. The 
momentum of relativistic electrons before scatter- 
ing is denoted by P. 

Since, according to the observational evidence, 
we observe the well localized pulses of high en- 
ergy emission, therefore the angle, must be very 
small. On the other hand, analyzing the excitation 
of oblique waves in a relativistic electron-positron 
plasma one can ar gue that the pitch ang le has to 
be extremely low ( Volokitin et a.1.1 Il985f ). Then, 
Eq. (fT2|) reduces to 

1 -0 



1 -f3 cos 6' + cos 6") 



(13) 



7mr 



We will study two principally different cases: (a) 
cos#' <C 1 and (b) cos 9' ~ 1. In the first case, we 
have 

(14) 



2 7 2 i + _^(i_cos0")' 

where we have taken into account the following 
approximate relation, 1 — (3 s» I/27 2 for (3 ~ 1. 
From this expression we see that for all physical 
quantities the frequency after the scattering is less 
than that of before the scattering and therefore, 
there is no possibility of increasing u> to the VHE 
band. 

By considering the second limit, Eq. (p~3j) re- 
duces to 

' (15) 



1 + ^(1 -costf") 



which, as in the previous case, leads to the similar 
result, lu' < ui. 

This investigation shows that the ICS cannot 
provide the VHE rad iation from the Crab pulsar 
detected by MAGIC (lAliu et alJl2008h . 



2.3. Curvature radiation 

Since particles are moving along the curved 
magnetic field lines continuously, they will emit 



the curvature radiation. On the other hand, in 
Eq. ([5]) we have neglected a term corresponding 
to the curvature emission. This in its turn, means 
that the following ratio 



V 



where ([Ruderman fc Sutherlandlll975l ) 



3h 3 c 

Ecur = TTlb X 5 
2 p 



(16) 



(17) 



and p is the curvature radius of magnetic field 
lines, must be less than one. By applying Eqs. 
(|llll6U7p . one can show that for typical magneto- 
spheric parameters of the Crab pulsar close to the 
light cylinder, the aforementioned ratio is negligi- 
ble only if the curvature radius exceeds the light 
cylinder radius, i?; c , approximately by three or- 
ders of magnitude. 

Generally speaking, since no contribution in 
emission comes from the closed magnetic field 
lines, we consider the open ones. On the other 
hand, in the dipolar field the region of almost 
straight field lines is just a tiny fraction of the 
emission area, leading to a negligible value of the 
high energy luminosity. 

If the beam component particles move along 
curved field lines, they experience the so-called 
curvature drift with the velocity: 



u b 



2^1 



(18) 



where LOu b — eBo/mc; and Bq is the background 
magnetic field induction. This velocity will even- 
tually create the drift current, J& r = en^Uf, which 
in turn, via the Maxwell equation 



47T 

(V X B) x = — J dr , 



(19) 



can create the toroidal magnetic field (by x we de- 
note a direction of the drift current. See Fig. [T}. 
This current is evidently less than the Goldreich- 
Julian (GJ) current Jgj — ent,c, since itb <C c. 
But on the other hand, the GJ current creates the 
corresponding magnetic field, B r rs AnJajRu/c, 
where R n « B{R)/B'{R) = R/3 (B' = dB/dR) 
is the length scale of the spatial inhomogeneity of 
the magnetic field. If we assume a dipolar configu- 
ration, then, by taking the value of the G J density, 
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Ufe Rj QB/(2wec), into account, one can show that 
the toroidal magnetic field equals jj^-B. Inside the 
light cylinder (R < R c ), this value is less than the 
background magnetic field- B and therefore such a 
toroidal magnetic field will be unable to rectify the 
twisted magnetic field lines. This implies that the 
curvature drift current, which is less than that of 
the GJ, cannot contribute to the process of recti- 
fying the field lines. However, in spite of that the 
drift current is not the source of the toroidal com- 
ponent, B r , it is a trigger mechanism for gener- 
atign the perturbed current, Ji = e(n®vl + riyUb) 
responsible for the creation of B r (see Eq. (|A10p ). 
where by upper script " 1" we denote the perturbed 
quantities. The source of the instability of current 
and the resulting magnetic field is the pulsar's ro- 
tational energy and the process is achieved via the 
parametrically excited curvature drift waves. The 
corresponding increment of the curvature drift in- 
stability can be presented by (see Appendix , for 
more details see (jOsmanov et al.ll2008l 120091) ): 



3 lo\ k x u b 
2 7fc k ec 



1/2 



J 



k x u b 

An 



kgc 



(20) 

where u>b is the beam component plasma frequency 
and 7bo - the Lorentz factor in an unperturbed 
state. k x and kg are the wave vector's components 
and £1 is the angular velocity of rotation. 

By considering the typical magnetospheric pa- 
rameters for the Crab pulsar close to the light 
cylinder, P sa 0.033s, j b Q = 10 8 and exam- 
i ning the perturbatio n lengthscale A w 10 8 cm 



( Osmanov et al. 2009), one can see that the in- 
crement is of the order of Is -1 . Comparing 
this value with the Crab pulsar's slowdown rate, 
4.2 x 10 _13 s _1 , we see that the instability growth 
rate exceeds by many orders of magnitude the 
slowdown rate, indicating that the mentioned in- 
stability is extremely efficient. 

It is worth noting that we have three types of 
the open field lines: (a) curved field lines which 
pass ahead of the rotation; (b) a tiny fraction of 
almost straight field lines and (c) curved field lines, 
lagging behind the rotation. 

If the initial perturbation of the toroidal mag- 
netic field satisfies the condition B r > 0, then 
such a perturbation will rectify all field lines which 
initially pass ahead of the rotation (suppose the 
clockwise rotation of the system) and will twist 



even more the magnetic field lines, which initially 
lag behind the rotation. In the case, B r < 0, the 
situation is opposite: the field lines initially lag- 
ging behind the rotation will be rectified. At this 
stage the curvature becomes infinity and as we see 
from Eq. (fl~8|) , the drift velocity tends to zero, 
saturating (killing) the instability. 

The investigation shows that, the curvature 
drift instability provides necessary conditions for 
an efficient mechanism of rectifying the field lines, 
leading to the negligible role of the curvature radi- 
ation, confirming our assumptions leading to Eq. 



3. Summary 

1. We have considered several emission mech- 
anisms for explaining the recently detected 
VHE emission from the Crab pulsar. 

2. Studying the synchrotron mechanism, we 
have shown that due to the cyclotron insta- 
bility efficiently developing on light cylinder 
scales, non-vanishing pitch angles are cre- 
ated, that leads to the efficient high energy 
synchrotron emission with the spectral index 
-1/2. 

3. The observational fact of the coincidence 
of signals in optical and high energy (> 
25GeV) intervals is in a good agreement 
with our model, in the framework of which, 
the cyclotron instability is excited in the op- 
tical spectra, which, via the synchrotron pro- 
cess leads to the high energy emission. 

4. Analyzing the inverse Compton scattering, 
we have found that for Crab pulsar's mag- 
netospheric parameters even very energetic 
electrons are unable to produce the observed 
photon energies. 

5. Considering the curvature radiation, we 
show that due to the curvature drift instabil- 
ity, the magnetic field lines are rectified very 
efficiently. This in turn, leads to a negligible 
role of the curvature emission process in the 
observed VHE emission along the aforemen- 
tioned rectified field lines. 

As we see, the detected coincidence of VHE 
and optical signals is an indirect confirmation of 
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the fact that (a) both spectra is produced by one 
source and (b) the only mechanism providing the 
detected high energy radiation is the synchrotron 
mechanism. This means that we observe the Crab 
pulsar towards these straight field lines, that is the 
reason why we do not see the curvature radiation 
coming from the twisted magnetic field lines. 
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A. Curvature drift instability 

In this section we study the process of rectifying the magnetic field lines due to the parametrically excited 
curvature drift instability. This instability is called parametric, because an external force - centrifugal force, 
plays a role of a parameter, changes in time and creates the instability. Generally speaking, the presence of 
an external varying parameter generates the plasma instability. The mechanism of energy pumping pr ocess 



from the external alternating electric fie ld into the electron-ion plasma is quite well investigated in (Silin 



1973 ; Galeev fc Sagdee"vlll973l ; lMaxlll973l ). In stead of considering th e altering electric field, one can examine 



the centrifugal force as a varying parameter ( Machabeli et al.l l2005). 



We start our consideration by supposing that the magnetic field lines are almost straight with very small 
nonzero curvature (see Fig. [l}. In this context we examine the field lines that are open, and thus have the 
curvature radius exceeding the light cylinder one, maximum by one order. Therefore, dynamics of particles, 
governing the overall picture of the curvature drift instability, can be studied, assuming that field lines are 
almost straight. In the framework of the paper we suppose that the plasma flow consists of two components: 
the plasma component composed of electrons and positrons (e ± ); and, the so-called, beam component (6) 
composed of relativistic electrons. It is well known that the dynamics of pla sma particles moving along the 



straight co-rotating magnetic field lines is described by the Euler equation: (Ma chabeli et al.ll2005f ) 



<>Va + (v a V)p Q - -c 2 7 a?V£ +— (e + -v„ x B ) , (Al) 



at m \ c 



where £ = ^/l — f2 2 i? 2 /c 2 , R is the coordinate along the straight field lines; p Q , v a , and e a are the momentum 
(normalized to the particle's mass), the velocity and the charge of electrons/positrons, respectively; a = 
{e ± , b} denotes the sort of particles and E and B are the electric and the magnetic field induction respectively. 
The continuity equation: 

y+V(« a v o ) = 0, (A2) 

and the induction equation: 

VxB = if + ^ £ J a , (A3) 
c at c *—i 

a—e^ 1 ,f> 

(where n a and Jq, are the density and the current, respectively) complete the set of equations for n, v, E 
and B. 

In the leading state the plasma is in the frozen-in condition: Eo + ^vo Q x Bo = 0, then, one can show that 
the soluti on to the Euler equation in Eq. (|A1[) for ultra relativistic particle velocities in the leading state is 



given by (jMachabeli k, Rogavalll994h : 

v Q e = w = c cos(Q,t + ip), (A4) 

where v., is the velocity component along the magnetic field lines and ip is the initial phase of each particle. 

For solving the set of Eqs. (|AHIA3[) . we wm linearize it assuming that, in the zeroth order of approximation, 
the flow has the longitudinal velocity satisfying Eq. (|A4[) and also drifts along the x-axis driven by the 
curvature of magnetic field lines (see Fig. [T]): 

u a = 2^L, (AS) 

UJ B a P 

where u a is the drift velocity; uJB a = ZaBo/mc; and Bq is the background magnetic induction. 
Let us expand the physical quantities up to the first order terms: 

+ (A6) 
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where "J = {n, v, p, E, B}. Then if we examine only the x components of Eqs. (|A1IA3|1 . and express the 
perturbed quantities as follows: 

tf^r) oc^f^exp^kr)] , (A7) 

by taking into account that kg <C k x and k r — 0, and bearing in mind that « cE^/Bq, one can show that 
Eqs. (|AlflA3|) reduce to the form: 

- i(k x u a + kgv^p 1 ^ = — v^Bl, (A8) 

^^-i{k x u a + k e v^nl l = ik x n a v l olx , (A9) 

- ikecBl = 4vr ^ e Q (n>^ + n l a u a ). (A10) 

a— e^ 1 .b 

According to the standard method (jOsmanov et al. I l2008h . after expressing v\ and in the following way: 

vl^v a y kA ^\ (ah) 

n^AT Q e ikA <*W, (A12) 

A ax (t) = ^( nt +^ + m sin[2 {nt + (A13) 

A ae (t) = ^sm(m + <p), (A14) 
and substituting them into Eqs. (|A8IA9|) . one can get the expressions: 

-ikA a (t') v ( t ')B T (t')dt' , (A15) 



mc7 Qo 



, = ie a n Q a k x ^ Aa{t) r l dt , r* e - ikAa (t") v ( t ")B r (t")dt", 
mc^ ao J J " 



e 



to Eq. (|A17p . the latter will simplify to the following form: 



a=e ± ,b cr=±l s,n,l,p 6 v ' 



(A16) 



which combined with Eq. (|A10[) , lead to the following form: 

-ik Q cBl(t)= V ^L e ikA °W [ e- ikA "^vJt')BJt')dt'+ 

a— ,d 

i £ JfL kx u a e lkA ^ f dt' f e- ikA "^v^{t")B r {t")dt", (A17) 



where w Q = e^J Aim Q a /m represents the plasma frequency. If we apply the following identity: 

,±ixsmy = J s (x)e ±isy , (A18) 



s 



xB r (u> + O (2 [a -l]+n-p + <j)) 



1 - 



k x U a 



LO - 



0(2s + n) 



E 



xe i<P(2[s-l]+n-p+a) , 

u 2 a k x u a ^ ^ J s (g a )Jn{h)Ji(g a )J P {h) 



E E 



iB r {uj + n (2[s -l + /j]+n-p+a))x e M^-i+^]+n- P +a) ^ 



(A19) 



where 



9a 



k g c 



4fi ' o 

and J s (x) (s = 0; ±1; ±2 . . .) is the Bessel function of integer order ( Abramovitz fc Stegan 19651 ). 

In order to solve Eq. (|A19[) . one has to examine similar equations, rewriting Eq. (|A19[) (with shifted 
arguments) for B r (uj ± O), B r (uj ± 20), etc.. This implies that we have to solve the system with the infinite 
number of equations, making the problem impossible to handle. Therefore, the only way is to consider the 
physics close to the r esonance condition, that p rovides the cutoff of the infinite row in Eq. (|A19[) ad makes 
the problem solvable (jSilin k, Tikhonchuklll970l ). 

Studying the resonance condition of Eq. (|A19|1 . one can derive the proper frequency for the CDI: 



w 



(A20) 



The present condition for physically meaningful case k x u a /2 < implies that 2s + n = and 2[s + /i] +n = 0. 

For solving Eq. (|A19[) . we examine the average value of B r with respect to p. Then, by taking into 
account the formula: 

^ J e iN *d<p = S N , , 

and preserving only the leading terms of Eq. (|A19|) . after neglecti ng the contribution f rom the plasma 
components, one can derive the dispersion relation for the instability (jOsmanov et al.l 120081 ): 



3u)%k x u b 



Jo 



40 



Jo 



kec 

7T 



2"/b kec 

By expressing the frequency as ui = ujq + iT it is easy to estimate the increment of the CDI: 

1/2 



3 w£ k x u b 
2 7bo k g c 



Jn 
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Fig. 1. — Here we show geometry in which we 
consider our system of equations. By e^, e r and 
e x unit vectors are denoted, note that e x _L e T ^. 
O is the center of rotation and C - the curvature 
center. 
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